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Abstract 

This work is the first of a series of articles addressing the energy optimization in ice 
hockey hahs. Here we outhne an analytic method to predict in which design and operating 
conditions the COP of the entire cooling system (refrigerator and cooling tower) COF sys 
is maximum. 

COF sys is investigated as a function of several variables, like electric consumption and 
brine physical properties. With this method, the best configuration and brine choices for 
the system can therefore be determined in advance. 

We estimate the optimal design of an average-sized ice rink, including pipe diameter, 
depth and brine type (ethylene glycol and ammonia). We also single out an optimal 
brine density and show the impact of the electric consumption of the pump on COF sys- 

Our theoretical predictions are validated with heat flow measurement data obtained 
at an ice hockey hall in Finland. They are also confronted with technical and cost-related 
constraints, and implemented by simulations with the program COMSOL Multiphysics. 

The multivariable approach here introduced is general, and can be applied to the 
rigorous preliminary study of diverse situations in building physics and in other areas of 
interest. 
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1. Introduction 

Indoor ice rinks are large buildings where an ice sheet (of area between 1000 and 
2300 m^) is cooled and kept at a temperature lower than -b^C by a complex refrigeration 
system, that includes compressors, brine pumps, evaporators and condensers. 

Clearly, the above implies a very high energy consumption, to which one should add 
also the heating of spectator stands and the air conditioning. A major challenge is then 
given by the implementation of methods for saving energy and reducing operational costs. 

This task can be accomplished in several ways. For instance, literature on mathemat- 
ical models for heat transfer in the concrete slab has recently flourished [ij [2j [S] H] [5] |6] , 
thanks also to the fast development of computational tools. 

However, the above mentioned research focuses on a precise study of heat transfer in 
the concrete slab, and/or on the ice surface. The results give a very precise analysis of 
the phenomenology of heat transfer, but do not manifestly provide with indications on 
how the ice hockey hall should be designed in order to save energy. 

Here we adopt a fairly different (yet complementary) approach, that aims to act to- 
ward this direction. The basic idea is the following: focus on the quantities of interest 
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and investigate their behavior in function of design parameterq^ This gives an optimal 
theoretical configuration, through a quick and simple mathematical analysis. After con- 
sidering technological constraints and costs, one eventually gives practical prescriptions 
on how the system should be improved (e.g. best brine, pipe size etc.). 

This prescription allows to single out the best design configuration of a system a 
priori^ still at the design stage. By means of a quick mathematical study, which reduces 
to a simple maxima/minima study of a multivariable function, one is able to predict the 
best design configuration without e.g. expensive empirical tests. 

This basic idea is so general that it can be used not only in any heat transfer processes, 
but also in any areas of interest to engineers. Moreover, as we show in this paper, 
it is certainly complementary to the accurate mathematical studies presented in e.g. 
[Il [HI m [To]. A detailed knowledge of the physical processes can and should be used 
indeed for enhancing the predictive power of the model. 

In this article we provide a specific example of this generic procedure, that is applied 
in this case to the Coefficient Of Performance (COP) of the entire cooling system (re- 
frigerator and cooling tower). This COPg^ys, here viewed as a multivariable function of 
the system parameters, such as pipe length and brine density, can be maximized by an 
appropriate choice of such parameters. 

A low efficiency of the cooling system (low COP) is a primary source of energy and 
money waste. Therefore the challenging topic of energy saving should include maximiza- 
tion of the COP. It is clear from the above that this reduces to simply studying the impact 
of each design parameter on the COP magnitude, thus finding the best configuration. 

We do not consider only the COP of the compressor because of the bound on the 
evaporation temperature (~ —15^(7) existing in ice hockey halls. The condensation 
temperature can be indeed changeqj according to the season and weather conditions. 
This cannot occur for the evaporation temperature, therefore there is not much freedom 
in the optimization of the compressor COP. We can instead circumvent this limitation 
by considering the full cooling system. This also constitutes a non-standard approach in 
this kind of analysis. 

Since the problem of energy saving in such a complex system as an ice hockey hall 
has several options, our study as a whole includes a number of perspectives divided into 
several articles. Here we address only the case of the COP as a multivariable function of 
the design parameters. 

The present paper is structured as follows. 

Section 2 discusses briefiy the experimental setup and the measurements. 

In Section 3, we consider the coefficient of performance COP^^s of the cooling system 
as a multi-variable function of the geometry of the ice rink, electric consumption and 
physical properties of the brine. We show an excellent consistency between predicted 
and measured COVgyg. 

In Section 4, the model proves that COV gys is much more dependent on the electric 
consumption of the cooling tower (related to condensation and evaporation temperatures) 
than on the pumping power. This agrees with known experimental results. 

In Section 5, we use our theoretical model to compare COV sys values according to 



-•^In this paper, the function is COV sys and the variables are pipe size, brine density etc. 

^This is done in the ice hockey hall in Leppavaara, where we use the condensation energy available. 
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Figure 1: The Leppavaara ice hockey hah (three compressors in the middle, "kentta 1" and "kentta 2" 
are the two ice tracks) 

different ethylene glycol concentrations and ammonia. We show that a more diluted 
ethylene glycol solution, or ammonia, are preferable. 

In Section 6 we focus on the design optimization, via simulations of heat transfer in 
the concrete slab with the program COMSOL. We show that the depth and size of the 
pipes inside the slab are not crucial to the heat transfer process. However, increasing the 
pipe number by 1/3 is an option which could be of interest in the analysis of control of 
the brine flow. 

In Section 7 the technological constraints of our theoretical predictions, based on the 
validated theoretical model, are discussed together with an estimation of costs. 

2. Ice hall structure and measurements 

In this paper, we consider a theoretical model based on measurements taken at an 
ice hockey hall built in 2009 in Leppavaara, Espoo, Finland (see Figjl]). It contains two 
ice tracks. 

Heat flux and temperature are measured with a plate installed at the ice-concrete 
slab interface. The sensors were placed about three meters from the rink so that during 
the resurfacing only one layer of warm water is spread over the heat flux plate. 

The heat flux plate and thermal sensors are connected to a logger, which records 
an input sent by the sensors every five minutes. Thus we get instantaneous values of 
the interface temperature in [C], and of the heat flux from the ice pad in [VK/m^], by 
multiplying the values saved in the logger by the individual calibration factor of the plate 
(see Figj2]). As we know the surface area of the ice track, we can calculate the total heat 
rate in [W]. 

In order to keep the ice frozen, the heat gains have to be canceled by the cooling 
power. This is generated with a refrigerating machine located in a separate container, 
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Figure 2: Heat flux and interface temperature in the ice track n.l in the Leppavaara ice hockey hall 



in a backyard of the main building. The electric consumption of the whole refrigerating 
machine including compressors, brine pumps, evaporators and condensers is measured 
by a pulse meter, which is connected to a logger set to save the cumulative value of the 
pulses received every five minutes (Figj3|. Again by multiplying this value we find the 
amount of used energy in [J], or dividing by this interval, the average power [W] during 
the period of measurement. 

The design power of the refrigerating machine is 610 kW. Refrigeration is guided by 
the thermal sensors of the ice track. The evaporation temperature is —15^(7 and the 
condensation temperature is between 24^(7 and 26^C There are three compressors, and 
the circulating refrigerant is NHs, otherwise ammonia. The brine pump for the ice track 
has electricity power 11 kW and the fiow is 28 1 /s. The brine used is an ethylene glycol 
solution at 35%. The power of the brine pump for the condensate is 7.5 kW. There are 
10 condenser^ and the power of the condenser fan is 3.2 kW. 



3. The COP as a multi-variable function 

Optimization of a cooling system can be achieved in several ways. Here we concentrate 
on maximizing the Coefficient Of Performance (COP) of the entire cooling system. It is 
defined as 

cooling power _ Q ... 
electric consumption Wp + i^/ 



^This class of cooling systems is very well described in [12] , 
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Figure 3: Total electric consumption Etot = Wp + E, as measured in four days. 



where we have spht the energy consumption in the two contributions of pumping power 
Wp of the brine pump and electric consumption E. The latter comprises the power 
consumed by the refrigeration system described in Section [2j namely condensers and 
compressors. 

Wp and the cooling power Q are instead functions of the volumetric flow V, brine 
temperature difference, geometry of the cooling system (length and diameter) and brine 
physical properties. 

Let us now expand the expressions for the cooling power Q and the pumping power 
Wp, to rewrite the COP in function of the variables of interest. First, expand the cooling 
power: 

Q mCpAT VAT 

Cpp 



COP, 



Wp^E Wp^E ""^Wp^E 



(2) 



Then, the total pressure drop is the sum of contributions from the collector pipes and 
from the pipes under the ice rink (subscript p). The pumping power can thus be written 
in the following form. 



Wp = VAp = V{f^^vl 



I £)24/5y4/5 



(3) 



where L and D are resp. the length and diameter of the main collector pipes, / and d are 
those of the 150 pipes under the ice rink and v^, Vj^^p are the mean bulk velocities of the 
brine in the collector and ice rink pipes, respectivel}Q / and fp are the corresponding 
friction coefficients. It can be shown that, for our average volumetric flow V = 37.57 //s. 



^Remember that Vm = V/A and Vm,p = Vp/Ap = V/{150Ap) 
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Figure 4: Ice pad construction at the Leppavaara ice rink. 



the fluid motion in the cohector pipes is turbulent (Re ~ 35000), while in the smaller 
pipes it is laminar (Re = 1052). This justifies the exponentials in Eqj3|. 

By substituting the above expression for the pumping power in (l2| and simplifying, 
we get our formula: 



COPgys — 



(4) 



Therefore one can express the COP implicitly as follows: 

COP = f{y,p,v,Cp,L,d,E), 



(5) 



namely by considering the COP as a function of the above variables, one could investigate 
the impact of each of them on the energy optimization. 

Finding for which values of these parameters the function COP is maximal corre- 
sponds to choosing the best possible configuration of brine fluid, volumetric flow, electric 
consumption and so on. We are therefore applying to our problem the method of func- 
tional optimization^ which is used in a number of diverse research topics [131 HU Ull HSl 

[171 [m [li Eo . 

In the next sections, we will take into consideration each and all of the variables 
entering Q, and calculate for which values the function ([5| is maximal. 
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Figure 5: Predicted COP, computed with Eq.(^, vs the COP as measured at the ice rink. 



4. Relative weight of pumping power and electric consumption in COPgys 

We first focus on the relative impact of Wp and E on tiie COP. We will see that 
this analysis is crucial for energy optimization, and that it leads also to the validation of 
Eq.Q with the experimental data. 

In Eq.([2|, the denominator Wp + E corresponds to the total electric consumption of 
the compressor, 

Etot = Wp^E, (6) 

that is measured at the ice rink in Leppavaara. This is, on the average, of order Etot ^ 
120kW (recall that the measurements were taken every 5 minutes). 

Since the data do not contain values of the pumping power and E separately, we 
compute V from the measured values of AT and from the cooling power Q. Substituting 
the according pumping power into Q and comparing the COP thus obtained with the 
measured value, we then find £^ as a fraction of the measured Etot- 

The above method returns the plot shown in Fig|5j showing an excellent agreement 
between data and theory and implying that E corresponds to 93.4% and Wp to 6.6% of 
the total electric consumption Etot- This is consistent with well-known results (see for 
instance [21 ). 

4-1- Impact of the volumetric flow on COF sys 

Consider Eq.Q as a function of only the volumetric fiow V and the temperature 
difference of the brine at inlet and outlet, AT. 
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Figure 6: COPsys as a function of volumetric flow and brine temperature difference at inlet and outlet. 

COP. ys is clearly proportional to AT, while the dependence on the volumetric flow 
is nontrivial. Thus we wonder whether there is a speciflc value of V such that the ratio 
cooling power vs pumping power reaches a maximum, leading to maximum coif] 

The COP can be plotted in a 2-dimensional graph as in Figure^ and as a function 
of only V if AT is flxed. We show this for AT = 1.5^C in Figure [7[^ where one can see 
that the COP increases until it reaches a maximum for a speciflc value of the volumetric 
flow, before decreasing again. 

The maximum corresponds to V ^ 901/ that is a very large value. However, this is 
expected as we have just demonstrated that the pumping power weights for only ~ 7% 
of the total electric consumption. In other words, we have shown that in principle, it 
is possible to reach the maximal COFsys by increasing the pumping power to very high 
values. 

Therefore we conclude that changing the magnitude of Wp, or in other words, working 
on the brine pump, is not determinant for energy optimization. This happens since 
technological limits to V do exist, and we have shown that the pumping power constitutes 
a small portion of the total electric consumption. 

Rather, looking at Eq.Q for the COP, our interest should be directed towards the 
following: 

1. Reduction of the electric consumption which comprises the chiller, condensers 
and heat exchanger. 

2. Control of the brine flow, in response to environmental factors on the ice rink, such 
as resurfacing and skaters. The aim is to minimize the (lowest) peaks in Fig js] and 
obtain a smoother shape of COFsys- 

3. Choosing the most efficient and cost-balanced brine for our conflguration. 



^Remember that E is independent of V. 
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as a function of only the volumetric flow. 



4. Verifying if it is possible to get an optimal design configuration (e.g. pipe diameter 
and number), to maximize the cooling efficiency. 

Points 1) and 2) require a dedicated analysis, therefore they will not be investigated 
here. We will address them in the forthcoming papers. In the next sections we focus 
instead on points 3) and 4). 

5. Brine choices 

First of all, recall Eq.Q, namely 



We want to see how the specific brine chosen is relevant for the COP. Let us study then 
how the above expression is dependent on the physical characteristics of the fluid. 

We have already shown that the pumping power counts for only ~ 7% of the total 
electric consumption, thus the kinematic viscosity u is irrelevant. One should then focus 
on the density p and the specific heat Cp. Since 



and in general Cp oc it is not trivial to decide which parameter should be increased 
to achieve higher COP. One then needs to write u and Cp in function of p, to obtain a 
plot COFsys vs p which singles out the point of maximum. 



COP, 



c^VAT 



sys — 




COPg^ys OC Cp , COPgys OC 



P_ 



(7) 
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Figure 8: Specific heat of an ethylene glycol Figure 9: Kinematic viscosity of an ethylene 
solution at 34% in function of density. glycol solution at 34% in function of density. 



However, there are no theoretical formulas u^Cp = /(p), therefore we interpolate 
experimental data taken at our laboratory, for an aqueous solution of ethylene glycol at 
34%. We obtain the following relations. 



u{p) = 24725p^ - 524575p + 278244 , 
Cp{p) = -2517.5p^ + 5337.6p- 2825.54. 



(8) 
(9) 



These provide with estimates which are sufficiently accurate for our purpose. The cor- 
responding experimental curves are reported in FigsjS] and [9j For typical values of our 
ice rink configuration {V = 37//s, AT - 1.5^C, E - 90kW), we then get the COP as a 



function of only the brine density as in Fig 18 

Last, let us study how COF sys changes according to the average brine temperature 
{Tin -^Tout)/'^' We use again our experimental data to derive the following interpolation 
formulas, 



p{T) = 1.1176 - 0.0002T, 
u{T) = O.OIT^ - 5.678T + 810.08, 
Cp{T) = -O.OOOIT^ + 0.0575T - 4.6061 . 



(10) 

(11) 
(12) 



To increase the precision of the fil[^ we consider the temperature interval of our interest 
263 < T < 278 K. The above equations are plotted in Fig|25l Fig(26] and Fig {27 
By plugging these expressions into Eq.Q, we obtain the COP values in function of the 
average brine temperature, for typical values of our ice rink configuration {V = 37 //s, 
1.5^C, Etot ^ 90 kW)^ as in Fig 21 This agrees with the experimental COP values. 



At 



5.1. Ethylene glycol concentration 

Let us go further with the study of the COP dependence on the brine characteristics 
according to Eq.Q. Here we aim to understand whether we can generalize the results 
of the previous section. To do this, we compare the values of COFgys obtained for four 
different ethylene glycol concentrations, with the same system configuration (volumetric 
fiow, electric consumption etc.). 



These fits are precise enough for our qualitative study 
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Figure 11: Simulation of heat transfer 
with three pipes. 



The results are plotted in Fig 21 it is clear that the lower the concentration, the 



higher the COP. To generalize this, we consider the plots in Figj22j Figj23] and Fig|24| 
which concern respectively the specific heat, density and kinematic viscosity of ethylene 
glycol in function of the concentration, at T = — lO^C 

By comparing these graphs with the COP in function of the concentration, we see that 
both V and p increase, while Cp decreases. However, GOV sys behaves like Cp. Accordingly, 
we conclude that, at least for an ethylene glycol brine, the specific heat is the most crucial 
physical parameter for optimizing the system COP. 



5.2. Ethylene glycol or ammonia? 

One last issue to take care of, is which brine would be better for maximizing COFgys- 
We consider an aqueous solution of 14% ammonia (NHs) and calculate the according 
COP with Q; the result is then compared to the COP with ethylene glycol {C2HQO2) 
at 20% and 34%. We find the following, 

COFsysiNHs) > COF sys{C2He02) , (13) 

for any ethylene glycol concentration. In particular, for a solution at 20% the gain is 9%, 
whilst for a solution at 34% the COP is larger by ~18%, that is quite remarkable. We 
conclude that according to our analysis, ammonia should be preferred to ethylene glycol. 



6. Design optimization 

We now address the question of how are the pipe number, pipe diameter and pipe 
depth in the slab relevant to the heat transfer. In other words: is it possible to obtain a 
more efficient ice cooling by altering any of these design parameters? 

In this case simulations are very handy, as they do not require expensive and time- 
consuming experiments, and provide already with good indications for further investiga- 
tions, if those are found to be necessary. 
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Let us consider the two-dimensional slab of Figj4j Following its symmetry, we focus 
on a 200 mm wide section, containing two pipes with 100 mm separation, and per- 
form simulations of conduction heat transfer with the program COMSOL Multiphysics 



(Fig 10). 



The setup is the following: the concrete slab has dimensions 200x120 mm, the pipe 
external diameter is 27 mm, the reference depth (at the pipe's edge) is 30 mm and 
the ice thickness is 30 mm as well. We used densities p = 2000 kg/ for concrete, 
p = 91S.9 kg/m^ for ice and p = 1.066 kg/m^ for a 34% ethylene glycol solution. 

A uniform heat flow q = 60VK/m^, consistent with our measurements, is applied to 
the ice surface, at initial temperature Ti = — 4^C Time-dependent heat transfer occurs 
between the ice and the pipes, at temperatures Tin = —9^C and Tout = — for the 
inlet and outlet pipe. These are the average tube temperatures in our data. The following 
results are obtained after t = 1200 s. 



6.1. Changing the pipe diameter 

First we consider different pipe sizes: d = 27 mm, d = 35 mm^ d = 40 mm. The latter 
diameter is certainly too large, but it is useful for visualizing more clearly the trend of 
the temperature vs the pipe size. 

Our results for the temperature at the ice surface and ice/concrete interface are 



plotted in Fig, 19 We find that, for an increase in the diameter from 27 to 40 mm, namely 
148%, the temperature changes by only ~ 1% (see Table [T]). Moreover, the temperature 
profiles do not change qualitatively, as shown by comparing Fig{T4]for d = 21mm and 
Fig{T2] for d = 40 mm. We conclude that the pipe size has a very small role in the 
optimization of the heat transfer. This is clearly due to the high heat capacity of ice and 
concrete. 



6.2. Changing the pipe depth 

The same method has been used when changing the pipe depth in the range 25 mm < 
h < 30 mm. Our results are reported in Table [l] and in Figj20] Also in this case, the 
system is not very responsive to a large increase of the parameter (here, the depth). 

From Table [T] we see indeed that pipes closer to ice by 20% lower the temperature by 
only 0.1 i^. Comparison of Fig 13 and Fig {l4] shows also that the temperature profiles 
are qualitatively the same. Therefore the pipe depth is basically uninfluent to cooling 
optimization. 

6.3. Two versus three pipes 



Depth (mm) 




Tcon{K) 


25 


267.97 


266.74 


27 


267.99 


266.78 


30 


268.03 


266.86 



Diameter (mm) 


T^ce{K) 


TconiK) 


27 


268.03 


266.S6 


35 


267.97 


266.70 


40 


267.91 


266.59 



Table 1: Maximum temperature of the ice surface and ice/concrete interface for different pipe depths 
and diameters. 



Consider the 200 mm wide slab in FigjToj with one inlet and one outlet pipe. We 
wonder how the heat transfer is modified by adding one more inlet pipe in the same 
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Width (mm) 



Width (mm) 



Figure 12: Temperature profile at the ice sur- 
face (solid) and ice-concrete interface (dashed), 
for two pipes with diameter d = 40 mm. 



Figure 13: Temperature profile at the ice sur- 
face (solid) and ice-concrete interface (dashed), 
for two pipes buried at depth h = 2b mm. 



width, thus instahing the pipes closer to each other. This is visuahzed in FigfTT] the 
distance between two contiguous pipes is now changed from 100 mm to 66 mm. 

COMSOL simulations show that the temperature profiles at the ice surface and 
ice/concrete interface are radically different in the two cases (compare also Figs 14 and 
[T5|. Besides the general lowering of T values at a specific the temperature profiles are 
much smoother and the heat transfer is more uniform for three pipes. 

Using 200 pipes instead of 150 for the whole ice rink could then be an interesting 
option. We see from Figs 14 and 15 a remarkable difference of T{x) at the ice surface. 
The temperature at x = 200 mm is almost 0.5^C smaller for three pipes, and its profile is 
much more uniform than in the case of two pipes. Moreover, with three pipes the whole 
system is clearly more responsive to temperature changes. This might be relevant when 
addressing control techniques of the brine flow, and will be taken into consideration in a 
further study. 
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Figure 14: Temperature profile at the ice sur- 
face (solid) and ice-concrete interface (dashed), 
in the case of two pipes. 



Figure 15: Temperature profile at the ice sur- 
face (solid) and ice-concrete interface (dashed), 
in the case of three pipes. 



7. Technological constraints and cost estimate 

The previous sections have provided with a series of insights and prescriptions for op- 
timization, which are based on the functional method we outhned (top-down approach). 
These are vahdated with measurements performed at an ice hockey arena. 

This section aims to complement such predictions with a bottom-up approach, dis- 
cussing their practical limits from the viewpoint of technical design and costs. We refer 
here to several sources, see for instance [22 | \23 \ [24 1 125], 

7.1. Pipe depth 

The refrigerating slabs contain pipes buried in concrete, therefore it is crucial to know 
the minimum depth of installation for the pipes under the upper surface of the concrete 
slab. 

There are several limitations on the thickness of the concrete layer. It must resist 
to a certain pressure (users, ice resurfacers), it works in compression and flexion. It 
should resist to fatigue, since several cycles of loading and unloading occur each day. 
The concrete must resist to multiple cycles of freezing and unfreezing, and must not let 
humidity reach the insulation layer (proper porosity), to avoid moisture-related problems. 

Some regulations set the minimum thickness of these types of slabs to 130mm, but 
these are for industrial applications (warehouses, sheds etc.). In our case, the maximum 
load is inferior to 10 kN/rn^ (distributed or concentrated) so the thickness of the slabs 
can be undersized. 

For an 80mm to 100mm polystyrene insulation layer (that can resist to maximum 
850 kN/m^), the minimum thickness of the concrete slab should be 65mm to resist to 
the maximum vertical load. Usually, this size goes up to 120mm in order to accommodate 
the pipes. Some North- American ice-rinks designs give us a thickness of Gin (150mm) 
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which is oversized. Moreover, at least one layer of wire mesh reinforcement is usually 
placed at the concrete-ice interface. 

Overall, the above limitations lead to design the concrete layer with a minimum 
concrete thickness of 25mm above the edges of the pipes. 

7. 2. Slab design and pipe size 

Different formulations of concrete can be used, like hemp concrete. A wire mesh 
reinforcement is needed especially if the arena is used for other purposes after the skating 
season is over, because the floor has to withstand a large spectrum of average load. In 
this case a strict regulation is to be applied to determine the distance between the steel 
wires and thus the layer thickness. One might also need to consider the heat transfer 
through this material. 

Note that when designing the concrete slabs, the stress and dilatations to which the 
pipes are submitted because of the variation of temperature should be taken into account 
as well. 

The PVC pipes are built to resist the pressure of the concrete layer above. In our case, 
a 27mm to 40mm diameter pipe will resist a pressure of a concrete layer with thickness 
from 25mm to 130mm (weight of the concrete plus steel reinforcements). 

Concerning the maximum depth of the pipes, the thickness of concrete under the 
pipes must not be too small to avoid direct contact with the insulation layer. However, 
as long as this contact is avoided there are no limitations. Placing the pipe in the middle 
of the concrete layer is a good solution (easy to set up). Note that another solution 
consists of prefabricating the slabs and seal them together on site, but this makes it 
more difficult to install a perfectly fiat surface. 

7.3. Price of the PVC Pipes 

The study performed in this paper is based on a 200mm wide slab containing 2 or 
3 pipes spaced out 100mm or 66mm each. Let us consider an Olympic ice rink of 30m 
width, for a total of 300 or 450 pipes. Then the total length of the PVC pipes would be 
17000m or 27000m (60m per pipe). Of course installing more pipes (3 per slabs rather 
than 2) will be more expensive, even if we chose the minimum diameter of pipes (the 
price- diameter relation is non-linear). What must also be taken into account is the price 
of the 180° bend connecting two parallel pipes. For N pipes we need N-1 bends. 

The average price for a 25mm to 40mm diameter PVC pipe for refrigerant applications 
is approximately 0.20€ to 0.50€per meter (industrial supplier). Prices on the market 
can differ from one supplier to another. Some PVC pipes are sold at 0.80€ per meter. 
The cost of a bend would be between 0.40€ to 0.80€. 

All in all, we find the following: 

• 2 pipes per slab: 17000m of pipes, 299 bends => 3 400€ to 15 000€ 

• 3 pipes per slab: 27000m of pipes, 449 bends => 5 600€ to 22 000€ 

It is clear that placing 3 pipes rather than 2 inside a 200mm slab clearly demands 
an additional expense of more than 50%. However, we will consider this option in the 
context of control in a future work, to determine whether the higher installation price is 
balanced by higher performance and lower maintenance expenses. 
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Figure 16: The COP as a function of the ethylene glycol concentration. 



7.4- Cost of the refrigerants: Ethylene Glycol vs. Ammonia 

The price of liquid Ammonia seems to be lower than the price of Ethylene Glycol, 
even if a serious increase of the costs of liquid Ammonia occurred some time ago. Today 
liquid Ammonia is sold indeed from $500 to $800 per ton (approx. 400€ to 600€), 
whereas Ethylene Glycol is sold at $800 to $1400 per ton (600€ to 1000€). 

Liquid Ammonia used to have the lowest price compared to all the other secondary 
refrigerants. For 1 kg, it was 5 to 8 times cheaper than others. It is also the lightest fluid 
used as a refrigerant, that makes it cheaper to use (from 11 to 17 times w.r.to. Ethylene 
Glycol). After the rise of prices, it is still 1.5 to 2 times less expensive (average price in 
April 2012: $500 per ton, ^380€). 

Moreover, even though the concentration of the refrigerant certainly determines the 
price, it is unlikely that in our case the use of Ethylene Glycol turns out to be more 
convenient than Ammonia. 

Overall liquid Ammonia seems to be the best solution (better COPsys^ many advan- 
tages such as chemical neutrality with pipes, olfactory detection of leaks, no green-house 
effect contribution and so on). 

7.5. Ethylene Glycol and Ammonia Concentrations 

According to the literature, excessive glycol concentrations should be avoided if we 
seek to maximize the heat transfer rate. Moreover glycol concentrations less than 20%vol 
shall not be used, since it is a nutriment source for bacteria at these low rates [26 . We also 
know that some existing cooling systems have the following features: brine at 34% vol 
Mono Ethylene Glycol, brine at 38%vol Mono Propylene Glycol, Liquid Ammonia at 
17-18%vol. Generally a cooling antifreeze liquid has about a 40%vol concentration. 

It is reported in Table [5] that for a -8°C freezing point (average pipe temperature in 
our data) the minimal concentration is 18,l%vol. This increases to 52%vol for -39° C [27] . 
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To sum up, ideal Ethylene Glycol concentrations should be in any case between 20 
and 52%vol. 

We are more likely to suggest using liquid Ammonia rather than Ethylene Glycol: 
another plea for this choice is that NH^ can be used with a lower concentration, which 
makes it even cheaper. 



Density (g/cm^) @20°C 


%vol 


Freezing Point (°C) 


1.024 


18.1 


-8 


1.030 


22 


-10 


1.067 


49.6 


-39 


1.072 


53.7 


-45 



Table 2: Ethylene Glycol concentrations and corresponding freezing points 



8. Conclusions and outlook 

In this paper we have outlined a functional method to study the optimization of 
cooling in ice hockey halls, and provided the optimal brine fluid and concentration for 
an average-sized ice rink in standard operating conditions. 

Investigating the COP of the entire cooling system as a function of several variables 
is an interesting option, as it provides with precise prescriptions on the system design 
that do not demand expensive and time-consuming preliminary tests. 

This a priori approach makes it possible to quickly address questions such as: the 
weight of volumetric flow and brine characteristics on the refrigerator performance; the 
role of the pumping power in the total electric consumption; choosing the most efficient 
and cost-balanced brine. 

We find that the pumping power Wp corresponds to ~ 7% of the total electric con- 
sumption, and that the volumetric flow is uninfluent on the system COPg^s for any 
technically acceptable values. 

Accordingly, we show that energy optimization should not be achieved by acting only 
on the pumping power, since a much more crucial role is given by the other sources of 
electric consumption such as the chiller, heat exchanger and condensers. Choosing the 
optimal brine fluid is also very relevant to achieve a high COFsys- 

By using the method of functional optimization, we have then singled out which brine 
fluid is optimal for the configuration at hand. It is shown that COF sys is mostly sensitive 
to the specific lieat of the fluid. Moreover, there is a specific value of the density for 



which the COP is maximal (Figure 18), regardless of the specific brine chosen. 

It follows that ammonia is a better choice than ethylene glycol, and that glycol 
at concentration between 20% and 34% is to be preferred. We have compared these 
predictions with technical contraints, and a cost estimation is provided in Section (T] 

Design possibilities regarding the pipe size and depth inside the concrete slab have 
been investigated as well. Transient heat transfer in a two-dimensional portion of the 
slab has been simulated with COMSOL Multiphysics, to give a consistent check to our 
assumptions at a qualitative level. 
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We confirm rigorously tiie well-known practical result that the pipe size and depth do 
not enhance appreciably the process. However, increasing the number of pipes provides 
with a more uniform temperature profile at the ice surface. This might be an interesting 
option to consider when assessing control methods of the brine flow. 

Our predictions are obtained with a model that is validated with measured data 
obtained at an ice rink in Finland, and are compared to technical and cost-related con- 
straints. 

We flnally stress that the multivariable approach here applied is general, and as such 
it can be applied to the rigorous study of diverse situations in engineering. 
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Figure 19: Temperatures for pipe diameters Figure 20: Temperatures for pipe depths 
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Figure 21: The COP as a function of the brine 

temperature for several ethylene glycol concen- Figure 22: Density of ethylene glycol as a func- 
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Figure 23: Specific heat of ethylene glycol as a Figure 24: Kinematic viscosity of ethylene gly- 
function of concentration. col as a function of concentration. 
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Figure 25: Density of an ethylene glycol solu- Figure 26: Kinematic viscosity of ethylene gly- 
tion at 34% in function of temperature. col at 34% in function of temperature. 




Figure 27: Specific heat of an ethylene glycol solu- 
tion at 34% in function of temperature. 
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